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Performance Analysis of Single- 
Phase Space Thermal Radiators 
and Optimization Through 

Taguchi-Neuro-Genetic Approach 


In the thermal management of spacecraft, space thermal radiators play a vital role as heat 
sinks. A serial radiator with proven advantages in ground applications is proposed and 
analyzed for space applications. From the performance analysis, specific heat rejection 
(SHR) of serial radiator is found to be higher than parallel radiator by 80% for 
maximum diameter of the tube, 47% for maximum thickness of the fin, and 75% for 
maximum pitch of the tubes under consideration. Also, serial radiator requires four times 
higher pumping power than parallel radiator with geometric parameters and a maximum 
mass flowrate under consideration. In serial radiators, the cross conduction between the 
fins has a significant effect on its thermal performance. Thus, conjugate heat transfer simu- 
lations and optimization operations are to be performed iteratively to optimize the serial 
radiator, which is computationally costly. To reduce the computational time, artificial 
neural network (ANN) is trained using conjugate heat transfer simulations data and com- 
bined with the genetic algorithm (GA) to perform optimization. Taguchi’s orthogonal 
arrays provided the partial fraction of conjugate heat transfer simulations set to train the 
ANN. Taguchi-Neuro-Genetic approach, a process that combines the features of three pow- 
erful techniques in different optimization phases, is used to optimize both parallel and serial 
radiators. The optimization aims to obtain a configuration that provides the lowest mass 
and lowest pumping power requirement for given heat rejection. Optimization results 
show that the conventional parallel radiator is about 20% heavier and requires about 
35% more pumping power than the proposed serial radiator. {|DOI: 10.1115/1.4052897] 


Keywords: conjugate heat transfer analysis, design of experiments, Taguchi’s orthogonal 
arrays, neuro-genetic optimization, multi-parametric optimization, extended surfaces/fins, 


heat exchangers, heat pipes, heat transfer enhancement 


1 Introduction 


Ever since USSR launched Sputnik in 1957, humans’ urge to 
explore space increasing day by day. With an increase in the 
number of human space flights, there is a need to develop new tech- 
nologies and optimize existing technologies. In a typical human 
space flight, the thermal control system (TCS), which is a subsys- 
tem of the Environmental Control and Life Support System 
(ECLSS), plays a vital role in maintaining the comfort cabin tem- 
perature for the crew and also in maintaining the electronic equip- 
ment within their operating temperature range. Heat sinks in TCS 
are usually the radiators, either body-mounted or deployed, which 
collect the heat from the coolant passing through the radiator 
tubes and reject it to deep space as electromagnetic radiation. 

Figure | schematically represents the simplified hydro-pneumatic 
circuit of the TCS. The air-liquid heat exchanger housed inside the 
cabin collects the heat dissipated in cabin air and transfers it to the 
primary coolant passing through TCS. Internal habitat loop (HL) 
consists of an air-liquid heat exchanger, electronics, interface heat 
exchanger, and the mechanically pumped fluid loop (MPFL), 
which includes pumps, solenoid valves, check valves, and flow pro- 
portioning valves. In the interface heat exchanger, heat is trans- 
ferred from the primary coolant to a secondary coolant of TCS. 
The secondary coolant is passed through the radiators using the 
MPEFL system to reject the heat to deep space, and this loop is 
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called the external radiator loop (ERL). This study aims at analyz- 
ing the performance of radiators of the TCS in ERL. 

From the early years of the space age to the recent past, research- 
ers have investigated the performance of heat pipe radiators and 
optimized the design parameters [1]. With increasing power 
requirements for the recent missions, better heat removal strategies 
with MPFL systems are required. MPFL systems, in general, are 
used for massive satellites with higher heat dissipation rates. 
Afshari et al. [2] performed the optimization of tube pitch and 
length for a fixed pipe diameter of MPFL parallel radiators using 
a genetic algorithm (GA). The literature on the radiator applications 
on the ground for providing thermal comfort to living spaces was 
reviewed. Salem et al. [3] studied the effect of Aluminum and 
Copper nanoparticles enriched fluid passing through the serial radi- 
ators for ground applications. Vosa et al. [4] experimentally ana- 
lyzed the heat rejection efficiency of parallel and serial radiator 
configurations. The authors concluded that the emission efficiency 
of serial radiator configuration is higher than that of parallel radiator 
configuration. Considering the higher emission efficiency of serial 
radiators, the present work proposes and studies serial radiators 
for space applications. 

Mass optimization of the MPFL parallel radiator was undertaken 
by Kumar Rai et al. [5] through a semi-analytical approach. Several 
optimization tools exist for designing radiators. Numerous authors 
reported using a GA, a powerful tool for optimizing nonlinear, dis- 
continuous, and nondifferentiable objective functions for optimiz- 
ing various thermal systems [6,7]. Dias and Milanez [8] brought 
out the capabilities of GA through a comparative study with work 
reported by da Silva et al. [9], where an extensive search optimiza- 
tion technique was used for the same problem of heat source 
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Fig. 1 Schematic of a hydro-pneumatic circuit of a thermal 
control system in a typical human space flight with parallel 
radiators 


location optimization on a vertical wall with natural convection heat 
transfer. The authors reported that the number of simulations per- 
formed by GA to obtain accurate results is only about 2% of 
those required in extensive search, thus, largely reducing the com- 
putational time. Thus, GA is selected as a tool for optimizing the 
parallel and serial radiator configurations. 

The thermal performance of parallel radiators can be mathemat- 
ically modeled through analytical expressions, making them amena- 
ble for direct optimization. Whereas for serial radiators, where the 
cross conduction between the fins plays a vital role, the thermal per- 
formance can only be accurately estimated through conjugate heat 
transfer simulations. Optimization of serial radiators requires the 
interconnection of conjugate heat transfer simulations with GA. It 
is an iterative process in which the genetic algorithm predicts the 
new population, for which conjugate heat transfer analysis is to 
be carried out to provide fitness values to each iteration population. 
Interconnection of conjugate heat transfer simulations with the 
genetic algorithm is computationally costly. The artificial neural 
network (ANN), which is trained using the conjugate heat transfer 
simulation data, estimates the output parameters for a given set of 
input variables. Leong et al. [10] brought out the capabilities of 
ANN to predict the thermal performance of plastic-ball grid array. 
The authors showed that ANN predicts the junction temperature 
with 21 times better accuracy and the thermal resistance with 69 
times better accuracy than multiple linear regression. Asok et al. 
[11] utilized ANN to relate the input parameters of static labyrinth 
seals with output parameters. The authors trained ANN using com- 
putational fluid dynamics data, and the trained ANN was used to 
predict the vortex loss coefficient and optimize the performance 
of square cavity labyrinth seals. Vortex generator for heat transfer 
enhancement in a plate-fin channel is optimized by Khan and Li 
[12] by combining the Bayesian-regularized artificial neural 
network, trained by computational fluid dynamics data and multi- 
objective optimization algorithm. Compared with the direct inter- 
connection of conjugate heat transfer simulations with GA, 
trained ANN and GA coupling is computationally faster. Combin- 
ing the capabilities of ANN and GA is termed as Neuro-Genetic 
optimization. Madadi and Balaji [13] optimized the location of mul- 
tiple discrete heat sources in a ventilated cavity using combined 
ANN and GA. ANN was trained using the data obtained from com- 
putational fluid dynamics, and the trained ANN was combined with 
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a genetic algorithm for performing the optimization. Lee et al. [14] 
used the results of finite element analysis of micro-compact heat 
exchangers for training ANN, and later ANN was combined with 
GA for performing optimization. Neuro-Genetic optimization is 
found to be efficient in terms of computational time and thus is 
selected for this study. 

The partial fraction set of conjugate heat transfer simulations to 
be carried out for training ANN is obtained from Taguchi’s orthog- 
onal arrays. The orthogonal array supports the design process and 
significantly reduces the number of simulations, saving time and 
resources. Jaya Krishna [15] optimized the phase change material 
heat sink using the Taguchi method, followed by melt fraction 
and maximum base temperature analysis. Aly [16] used the 
Taguchi method for finding the optimum condition for the heat 
transfer design parameters of the coiled tube-in-tube heat exchanger 
under turbulent flow conditions. To observe the effect of specific 
geometric parameters of H-type finned tube heat exchangers, 
Wang et al. [17] used Taguchi’s process. The effects of Reynold 
number, fin height, and pitch on the heat transfer efficiency of a 
heat exchanger having cylindrical pin fins positioned in a rectangu- 
lar channel were experimentally studied using the Taguchi method 
by Sahin [18]. Altan [19] used the L>7 orthogonal array to produce 
21 training data and six test data for the shrinkage prediction algo- 
rithm of a back-propagation neural network. Taguchi orthogonal 
arrays possess considerable advantages in providing the partial frac- 
tion set of conjugate heat transfer analysis required to train ANN 
and are adopted in this study. 

The present paper combined the advantages of three efficient 
methods in various phases of the optimization process, namely, 
Taguchi method for the design of experiments, the trained artificial 
neural network using conjugate heat transfer simulations, for esti- 
mating output parameters for a set of input variables, and the 
genetic algorithm for carrying out multi-parametric optimization. 
Combining the Taguchi method, ANN, and the GA for optimization 
studies of complicated problems is called Taguchi-Neuro-Genetic 
methodology. 

Summarizing this section of the paper, the literature survey on 
space radiators indicates that a large portion of the literature avail- 
able in the public domain deals with optimizing parallel radiator 
configuration for space applications. The motivation behind the 
present work is to analyze the performance of the serial radiator, 
which has proven advantages on the ground, for space applications 
and to arrive at an optimized serial radiator configuration. The spe- 
cific objectives of the study are as follows: 


(a) To propose and study the serial radiator configuration for 
space application. 

(b) To analyze the effects of various geometric and flow param- 
eters on the performance of an MPFL parallel/serial radiator 
through conjugate heat transfer simulations. 

(c) To optimize the parallel and serial radiator configurations 
using Taguchi-Neuro-Genetic approach 

(d) To bring out the superior radiator configuration for different 
optimization objectives 

(1) Case 1: the objective is to minimize the sum of nondi- 
mensional mass and pumping power for a given heat 
rejection 

(11) Case 2: the objective is to minimize nondimensional 
mass for a given heat rejection 

(111) Case 3: the objective is to minimize nondimensional 
pumping power for a given heat rejection. 


2 Radiator Configurations and Inputs for the Analysis 


2.1 Radiator Configurations. For the present study, radiators 
are envisaged as four flat panels arranged around the service module 
of a human space flight, as shown in Fig. |. Two radiator configu- 
rations are considered for the study: (1) conventional parallel radi- 
ator and (2) proposed serial radiator. 
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Fig. 2 Parallel radiator configuration 
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Fig. 3 Serial radiator panel configuration: (a) schematic and (b) 
realized radiator 


2.1.1 Parallel Radiator Configuration. In the parallel radiator 
configuration, the tubes inside the radiator are arranged parallel to 
each other, as shown in Fig. 2. Flow is drawn to the tubes from a 
common inlet header and dumped into a common outlet header. 
As the radiators are segregated as four parallel independent 
panels, the flowrate to each radiator panel is 1/4th of the total flow- 
rate of the external radiator loop of the thermal control system. With 
a pitch of 100 mm, a radiator panel of width 1 m consists of ten 
tubes in parallel. Thus, the coolant flowrate to each tube for this 
case is 1/40th of the total coolant flowrate. 


2.1.2 Serial Radiator Configuration. In the serial radiator con- 
figuration, the tubes in the radiator panel are connected in serial. 
Thus, each radiator panel is considered to have one tube running 
all along the panel in a serpentine fashion, as shown in Fig. 3. 
The coolant flowrate to this tube is 1/4th of the total coolant flowrate 
of the ERL of TCS. Unlike parallel radiator configuration, serial 
radiator configuration does not require common inlet and outlet 
headers. 


2.2 Parameters Influencing the Performance of a 
Radiator. The performance of the radiator is assessed through 
the following parameters: 


(a) Heat rejection to space (Q) in Watts is given by Eq. (1) 


QO = o&FA(To4 — Tonk) (1) 


Journal of Thermal Science and Engineering Applications 


(b) Total mass of the radiator (/) in kg is given by Eq. (2) 
M=m, +m +m3+m4+mMs5 (2) 
(c) Specific heat rejection (SHR) in W/kg is given by Eq. (3) 
SHR =Q/M (3) 
(d) Pumping power (II) in Watts is given by Eq. (4) 


Total Volume flow rate x Total pressure drop 





Efficiency of pump 


1 AX aor/ Pp AP + 15,000) 
"p 


(4) 


Since the pump has to circulate the coolant in the entire ERL and 
as four radiators are connected in parallel in ERL, the total mass 
flowrate in ERL is four times the mg; per radiator. Apart from pres- 
sure loss in the radiator (AP), the pump in the ERL has to overcome 
the frictional losses, bend losses in the main plumb lines as well as 
pressure drop in components in the ERL such as interface heat 
exchanger, solenoid valves, check valves, and flow proportioning 
valve. The losses in main plumb lines and component losses are 
independent of the geometric configuration of the radiator tubes 
and thus considered as a constant value of 15,000 Pa. Thus, the 
total pressure loss that the pump has to overcome is (AP+ 
15,000) Pa. The overall efficiency of the pump (np) is considered 
as 30%. 

For an ideal radiator, the heat rejection requirement shall be met 
with the minimum mass of the radiator and the pumping power 
requirement. 

For a given external environment and radiator area, the parame- 
ters influencing the performance of the radiator, as depicted in 
Fig. 4, are as follows: 


(1) Geometric parameters 
(a) Diameter of the tube running through the radiator panel 
(di) 
(b) Thickness of radiator fins (¢) 
(c) Pitch of the tubes in radiator panel (P;) 
(d) Thermo-optic properties of the radiating surface 
(ii) Flow parameters 
(e) Fluid type and its properties 
(f) Mass flowrate of fluid passing through the tubes (mgr) 
(¢) Temperature of fluid at the inlet of the radiator (T;,) 
(ii) Fin and tube material properties 


For the present study, the thermo-optic surface properties of the 
radiator, fluid type, fluid inlet temperature to the radiator, fin and 
tube materials are frozen, keeping d;, ty, P;, and mao; aS variables. 


2.3. Material Selection. Aluminum alloy AA6061-T6 is 
selected as the fin and tube material for better weldability, brazeabil- 
ity, workability, and machinability. Strength and corrosion resis- 
tance are increased with tempering. 

3M NOVEC ~ 7100 is selected as the heat transfer fluid for its 
better thermophysical properties and lower freezing point. Table | 
gives the thermophysical properties of the selected materials. 





Fig. 4 Schematic of radiator showing variable parameters 
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Table 1 Thermophysical properties of selected materials 


3M NOVEC™ 
AA6061-T6 [20] 7100 [21] 
Density (p), kg/m° 2700 1510 
Conductivity (K), W/mK 160 0.069 
Specific heat (C,,), J/kgK 896 1183 
Dynamic viscosity (uw), Pas — 5.8x10* 
Freezing point, °C — —135 


2.4 Inputs for the Analysis:. This study considered an orbit of 
about 400 km altitude and with an attitude of spacecraft such that 
there is no direct solar flux on the radiators. In such orbit, the 
albedo (A) is 239 W/m? and earthshine (E) is 200 W/m? [22]. 

Typically, the radiators are coated with white paint to increase the 
radiation emission to deep space and decrease external heat load 
absorption. Typical surface optical properties of white paint are 
solar absorptivity of 0.15 and IR emissivity of 0.9. Thus, the exter- 
nal heat load on radiators is given by Eq. (5) 


dext =A; X A+epR XE 
Gext = 0.15 x 239 + 0.9 x 200 


”. dext = 216 W/m? (5) 


While on the radiator toward the space side, the external heat load 
remains zero. Thus, an average external heat load on all four radia- 
tors will be in the order of 100 W/m” and this value is considered 
for the present study. 

The external heat load and the ambient space temperature are 
converted to equivalent sink temperature through Eq. (6) [23] 


1 
Tsink = T space + (ext/ F oe) 4 (6) 


3. Methodology 


The methodology adopted for optimization of radiators is 
Taguchi-Neuro-Genetic approach, which is depicted through a 
block diagram given in Fig. 5. The radiator performance is analyzed 
through the steady-state conjugate heat transfer simulations. The 
minimum partial fraction set of simulations to uniformly cover 
the entire design domain is obtained through Taguchi’s orthogonal 
arrays. Conjugate heat transfer simulations are then carried out for 
the set obtained through Taguchi’s orthogonal array. The data 


Artificial Neural Network © 
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obtained from the conjugate heat transfer simulations is used to 
train the ANN, and the trained ANN is coupled with GA to optimize 
radiators. Detailed briefing of each process in the methodology is 
given below. 


3.1 Conjugate Heat Transfer Simulations. The conjugate 
heat transfer simulations are carried out using the “Electronics 
System Cooling” module of Siemens’ NX software. A fully 
coupled pressure—velocity solver with a mixing length turbulence 
model is used for the simulations. The fluid is modeled as 3D 
HEX elements with boundary layer mesh at the fluid—solid inter- 
face, and solid is modeled as 3D TET4 elements. Fluid mesh param- 
eters are given as 0.5 mm element size with a surface maximum 
growth rate of 1.3 to obtain a wall y+ of 0.9. Figure 6 displays 
the mesh used for conjugate heat transfer simulations. The fluid 
and fin elements are provided with an advanced coupling option 
to ensure the convective heat transfer from fluid elements to neigh- 
boring fin elements. 

In parallel radiator configuration, the coolant to all parallel tubes 
enters at the top and exits from the bottom. Thus, the significant 
temperature gradient on the radiator surface is longitudinal (top to 
bottom). In the lateral direction, about the centerline between two 
tubes, the temperature distribution is symmetrical. Hence, an adia- 
batic boundary condition can be considered at the centerline 
between the adjacent tubes. Thus, it is possible to treat the 
problem using the fact that the analysis of one set of finned-tubes 
can be extended for the whole radiator. 

While in serial radiator configuration, the exit of one straight tube 
is connected to the inlet of the other straight tube through U-bent 
tubes. The temperature distribution about the centerline between 
the tubes in the lateral direction is not symmetrical. Thus, for 
serial radiator configuration, the entire panel has to be modeled 
for the simulation. 


3.1.1 Assumptions 


(1) Fin material is assumed to be isotropic. 

(ii) The thermophysical properties of fin material and coolant 
are assumed to be constant within the temperature range 
of the problem. 

(111) Orbit average external heat load is used for estimating the 
average sink temperature. 


3.1.2 Boundary Conditions 


(i) Radiation from the flat surface of the radiator to sink 
temperature. 

(11) Adiabatic boundary condition is specified on all other exter- 
nal surfaces of the radiator. 


Genetic Algorithm 
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Fig. 5 Block diagram of Taguchi-Neuro-Genetic Approach 
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Fig.6 Mesh with fluid and solid elements of radiator (gray color 
is fluid) 
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Fig. 7 Configuration chosen for validation studies 


(111) Coolant mass flowrate and inlet temperature are specified at 
the inlet of the tube. 

(iv) Pressure outlet is specified at the outlet of the tube. 

(v) Thermal coupling between fluid and fin elements. 


3.1.3 Validation of Conjugate Heat Transfer Simulations. The 
baseline configuration for validating the conjugate heat transfer 
analysis is chosen from the work carried out by Kumar Rai et al. 
[5S]. A radiator with 20 numbers parallel finned-tubes shown in 
Fig. 7 is modeled as per the dimensions given in Table 2. Kumar 
Rai et al. [5] provided the properties of materials to be used for com- 
parative studies. External absorbed heat load (50 W/m?) on the 


Table 2 
et al. [5] 


Input parameters of the validation case by Kumar Rai 


Geometric parameters 


Fin thickness, t, 0.6 mm 
Inner diameter of the tube, d; 7.5mm 
Outer diameter of the tube, d, 8.5 mm 
Length of fin, Lp 48.4 mm 
Length of tube, Lr 475 mm 
Thickness of tube, t, 0.5 mm 
Number of tubes, N, 20 
Input conditions 

Mass flowrate of coolant, mgt 0.002 kg/s 
Fluid inlet temperature, 7;, 325 K 
Emissivity of radiating surface, € 0.8 
Space ambient temperature, T pace OK 
External absorbed heat load on top and bottom surfaces, diop 50 W/m” 
and bot 

Sink temperature on top and bottom surfaces, Toink 182K 
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Table 3 Results of validation studies 


Validation case by 
Kumar Rai et al. [5] 


Conjugate heat 
transfer analysis 


Fluid outlet temperature, K 317.09 317.65 
Mass of radiator, kg 2.265 2.2645 
Heat rejection capacity, W/m? 588.61 588.11 


radiator is converted to sink temperature (182 K) as per the proce- 
dure provided in Sec. 2.4. Mass flowrate and inlet temperature of 
the fluid, as provided in Table 2, are chosen as input parameters. 
Table 3 compares the results of conjugate heat transfer simula- 
tions and the results reported by Kumar Rai et al. [5]. The surface 
plot of temperature is shown in Fig. 8. It is observed that the differ- 
ence is less than 1%, and thus, the present methodology is validated. 


3.2 Taguchi Orthogonal Arrays for Design of Experiments/ 
Simulations. Ideally, a full factorial design will identify all possi- 
ble combinations for a given set of factors. With number of vari- 
ables as 4 and each variable has a level of variance of 5, the full 
factorial simulations would require 5* = 625 runs. Conducting full 
factorial of 625 runs is computationally time consuming. To 
reduce the number of simulations to a practical level, Taguchi’s 
orthogonal arrays for designing experiments (in this study these 
are simulations) are utilized. The columns of designed orthogonal 
arrays are balanced in two ways. First, the columns are balanced 
so that they have equal levels for each factor, 1.e., the columns 
are balanced within themselves. Second, any two columns are 
balanced among themselves such that together an equal number 
of possible combinations are generated. Design of experiments 
involves the following steps: 


(a) Selection of independent variables: The number of variables 
(n) in this study is 4; inner diameter and pitch of tubes in the 
radiator, thickness of radiator fin, and mass flowrate of 
coolant in the radiator 

(b) Selection of the number of level settings of each variable: 
Each variable is varied at a level (1) of 5 to capture the non- 
linear trend of the output parameters. 

(c) Selection of orthogonal array: For selecting the orthogonal 
array, the minimum number of experiments to be conducted 
shall be fixed based on the degrees of freedom, which is 
given by Eq. (7) 


Nragucri =1+ YG -1) (7) 

i=1 
From Eq. (7), Nraguchi 18 17 for n=4 and /=5. The orthogonal 
array shall be greater than the degrees of freedom and shall be 
second power of level of variance. Orthogonal array L25, which 
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Fig. 8 Surface temperature plot for validation case 
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Table 4 


Level 1 
Diameter of the tube (d;), mm 4 
Thickness of fin (t), mm 0.6 
Pitch of tubes (P;), mm 20 
Mass flowrate per radiator (770), g/S 10 


meets both the requirements is the right choice for this study. 
Table 4 provides the input variables and their levels and Table 5 
provides the Lz5 orthogonal array. Columns of the Table 5 represent 
variables and the number in the orthogonal cell represents the level 
assigned to the variable. Each row represents the combination that 
serves as the simulation set. For instance, simulation no. 10 corre- 
sponds to d;=8 mm, t= 1.6 mm, P;= 20 mm, and mg = 25 g/s. 


3.3. Artificial Neural Network to Estimate Output 
Parameters for a Set of Input Variables. Artificial neural 
network is created through a standard Neural Net Fitting toolbox 
available in MATLAB version R2020a. The conjugate heat transfer 
analysis data set is used for training the ANN. The ANN with 
input and output layers and one hidden layer, as shown in Fig. 9, 
is used in the present investigation. The number of neurons in the 
hidden layer is a tradeoff between the flexibility and overfitting of 
the network. Higher the number of neurons better will be the flex- 
ibility. With higher number of neurons, the total parameters (includ- 
ing weights, bias, etc.) in the network will be higher. Because of the 
fact that the total training data set shall be larger than the total 
parameters to avoid overfitting, the number of neurons shall be 
kept minimum. For the present study, where four input parameters 
and three output parameters exist, the number of neurons in hidden 
layer is selected as 4. The total number of parameters is 35 (Weights 
from the input layer to hidden layer: 4 x 4 + bias on hidden neurons: 
4+ weights from hidden layer to output layer: 4 x 3 + bias on output 
parameters: 3). Thus, training data shall be greater than 35 cases, 
whereas the minimum set of simulations obtained from Taguchi’s 
orthogonal array is only 25. Additional random simulation sets 
are conducted using conjugate heat transfer simulations and are 
used to train the ANN. In addition to increasing the number of simu- 
lations, the Bayesian regularization algorithm is used to avoid over- 
fitting and increase the generalization of the network. All the input 
variables and output parameters are normalized between —1 and 1, 
and the hyperbolic tangent activation function is used at each 
neuron. About 85% of the conjugate heat transfer simulation data 
is used for training, and 15% is used for test. 

The regression obtained for the network is close to 0.995, with a 
mean squared error of less than 0.001, indicating that the trained 
ANN predicts output very close to the target value for given input 
variables; the same can be observed from the plot given in 
Fig. 10, where the comparison is drawn between the heat rejected 
estimated by conjugate heat transfer analysis and that predicted 
by ANN. Trained ANN predicts the heat rejection within 1% of 
the value estimated by conjugate heat transfer analysis. 


3.4 Genetic Algorithm for Multi-Parametric 
Optimization. The trained ANN is coupled with the GA to 
perform the multi-parametric optimization. Standard GA function 
in MATLAB® is used. The initial population of random chromosomes 
containing the input parameters, d;, t, P;, and mg, generated by 
GA, is fed to ANN to predict output parameters; total heat rejected, 
mass, and pumping power required. Based on the objective func- 
tion, the fitness score is allocated to each chromosome. A cross-over 
operation between highly fit and inferior chromosomes 1s carried 
out, followed by a mutation operation, in which the genes of one 
of the individual chromosomes are mutated to form the new gener- 
ation population. The new generation population is then fed to the 
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Input variables and their levels 


Level 2 Level 3 Level 4 Level 5 
8 12 16 20 
0.8 2 1.4 1.6 

60 100 140 180 
25 50 75 100 


trained ANN to obtain the objective function based on the out 
parameters. The iterative loop continues until the output parameters 
converge. The flowchart for the Taguchi-Neuro-Genetic optimiza- 
tion process is displayed in Fig. 11. 

The optimization objective is to obtain a combination of input 
parameters that provide the lowest mass and lowest pumping 
power required for a given total heat rejection. For achieving this, 
a composite objective function, which is the sum of nondimensional 
mass and nondimensional pumping power, is used 


Composite Objective function = (M/M,) X w, + (II/II,) X w2 (8) 


4 Results and Discussion 


The performance of parallel and serial radiator configurations 1s 
independently assessed through the conjugate heat transfer simula- 
tions and compared. The baseline configurations for both parallel 
and serial configurations, as shown in Table 6 are selected to 
carry out the performance analysis. Apart from the simulation sets 
obtained through Taguchi’s orthogonal arrays, nineteen other simu- 
lations were carried out, for each of parallel and serial radiators, by 
varying one of the parameters for each simulation with respect to 
the baseline configuration. 

Figure 12 displays the comparison of parallel and serial radiator 
performance obtained by varying the tube inner diameter. For both 
the configurations, heat rejected by radiator decreases with an 
increase in tube diameter. This is because, for a given mass flowrate, 
with an increase in tube diameter, the velocity of flow in tube 
decreases, resulting in lower Reynolds number and eventually 
lower heat transfer coefficient from fluid to the fin surface. It is 
also observed that, for parallel radiator configuration, the increase 
in tube diameter beyond dj/d;, of 0.4 has an insignificant effect 
on the heat rejection because of the fact that at dj/d;, of 0.4 and 
beyond, the flow in the tube remains in the laminar regime. Rey- 
nolds number for a tube diameter dj/d;, of 0.2 1s 2745 (transition 
regime), while for a tube diameter dj/d;, of 0.4 is 1373 (laminar 
regime). It is also evident that the mass of the radiator increases 
with the increase in tube diameter for both the configurations. 
The findings on the variation of heat rejection and the radiator 
mass with tube diameter agree with the study reported by Afshari 
et al. [2]. The increase in the mass of parallel radiator is steeper 
compared with that of the serial radiator due to the increase in 
header mass of parallel radiator at a higher tube diameter. 
Figure 13 displays the comparison of SHR of parallel and serial 
radiators with variation in tube diameter. SHR is defined as the 
ratio of total heat rejection to the total mass of the radiator. A 
decrease in heat rejection and increase in mass leads to a decrease 
in SHR with increased tube diameter for both parallel and serial 
radiators. It can be noted that, at any tube diameter, the SHR of 
the serial radiator is higher than the parallel radiator indicating the 
serial radiator is a better thermal performer than the parallel radiator. 

Figure 14 displays the comparison of parallel and serial radiator 
performance with varying fin thickness (¢). With an increase in fin 
thickness, a gradual increase in heat rejection is observed for the 
both configurations. The cross-sectional area available for heat con- 
duction increases with fin thickness, leading to a higher average fin 
temperature. The fin tip temperature will be closer to the fin base 
temperature for higher fin thickness, leading to higher radiator effi- 
ciency. The mass of both the radiator configurations increases 
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Table 5 Standard Lz, orthogonal array 


Mass flowrate 


Simulation Diameter of Thickness of Pitch of per radiator, 
no. tube, d; fin, ty tubes, P; Mot 
1 1 1 1 1 
2 i 2 2 2 
3 1 e. 3 3 
4 1 4 4 4 
5 1 5 5 5 
6 2 1 2 3 
7 2 2 3 4 
8 2 5 4 2) 
9 2 4 =) 1 
10 2 | 1 2 
11 3 l 3 i) 
12 3 2 4 1 
13 5 3 5 2 
14 3 4 1 3 
15 3 5 2 4 
16 4 l 4 2 
17 4 2 5 2 
18 4 3 1 4 
19 4 4 2 5 
20 4 5 3 1 
21 5 l 5 4 
22 5 2 1 5 
23 5 3 2 1 
24 5 4 3 2 
25 5 5 4 3 


linearly with an increase in fin thickness. Figure 15 displays the 
comparison of SHR of parallel and serial radiator with variation in 
fin thickness. The increase in mass dominates the increase in heat 
rejection, leading to a decrease in SHR with an increase in fin thick- 
ness for both parallel and serial radiators. It can be noted that the 
SHR of the serial radiator is higher than the parallel radiator at 
any fin thickness. 

Figure 16 displays the comparison of parallel and serial radiator 
performance with varying pitch of tubes (P;). With an increase in 
pitch, the heat rejection by both the configurations reduces. At the 
pitch value (P;/P;,) of 0.11, the heat rejected by the parallel and 
serial radiator configurations is of the same order because, at very 
low pitch, the significant portion of the radiator area is covered 
with the tube base resulting in better heat rejection to space and 
thereby, the insignificant difference exists between heat rejection 
by parallel and serial radiators. The lower the pitch, the closer is 
the fin tip temperature to that of the fin base. Hence, low pitch 
leads to higher heat rejection capacity and higher efficiency of the 
radiator. As the pitch of the tubes increases, the reduction in heat 
rejected for the parallel radiator is steeper than that of the serial 
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Fig. 10 Comparison of results between conjugate heat transfer 
analysis and trained ANN 


radiator. Mass of both the radiator configurations reduces with an 
increase in the pitch of tubes. Figure 17 displays the comparison 
of SHR of parallel and serial radiators with variation in pitch of 
the tubes. At smaller pitch of the tubes, up to P;/P;, of 0.55, with 
increase in pitch of tubes, the decrease in mass is prominent than 
the decrease in heat rejection, leading to a higher SHR for both par- 
allel and serial radiators. At higher pitch of the tubes, above P;/P; , 
of 0.77, the decrease in heat rejection is compensated by the 
decrease in mass, leaving SHR almost unchanged with the increase 
in pitch of tubes. At any pitch of tubes, it can be noted that the SHR 
of serial radiator is higher than the parallel radiator. 

A comparison between parallel and serial radiator performance 
with change in the coolant mass flowrate is displayed in Fig. 18. A 
higher mass flowrate leads to a higher velocity of flow inside the 
tube, resulting in a higher Reynolds number and higher heat transfer 
between the flow and the fin surface. It is observed that the profile of 
the heat rejection curve for parallel radiator changes at a mass flow- 
rate (Mgot/Maot,o) Of 0.75, because the flow inside the tube transits to 
turbulent regime above (Mqo1/Maot,o) Of 0.75. The increase in the pres- 
sure drop with increase in the mass flowrate results in a higher 
pumping power requirement for both configurations. Flow transition 
from laminar to turbulent regime occurs for the parallel radiator con- 
figuration, whereas the flow for the entire range of flowrates of the 
present investigation remains in the turbulent regime for serial radi- 
ator configuration. For the coolant mass flowrates under consider- 
ation, the heat rejected and pumping power required by the serial 
radiator are higher than those of parallel radiator. 

Serial radiator provides higher heat rejection at lower mass but 
requires higher pumping power for a given geometric and flow 
parameters, than parallel radiator configuration. However, a detailed 
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Fig. 11 Flowchart of optimization through Taguchi-Neuro- 
Genetic approach 


Table 6 Baseline configuration used for parallel and serials 
radiator configurations 


Parallel Serial 
Geometric parameters 
Tube internal diameter (d;), mm 8.0 8.0 
Length of the radiator (L), m 1.5 1.5 
Width of the radiator, m 1.0 1.0 
Fin thickness (¢), mm 1.2 12 
Pitch of tubes, mm 100.0 100.0 
Flow parameters 
Mass flowrate of the fluid in the radiator (mot), g/s 50.0 50.0 
Temperature of coolant at the inlet of the tube (7;,), K 283 283 
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P-M: Mass of parallel radiator; 
S-M: Mass of serial radiator; 


P-Q: Heat Rejected by parallel radiator; 
S-Q: Heat Rejected by serial radiator; 


Fig. 12 Comparison of performance of parallel and serial radia- 
tors with variation in tube diameter 
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Fig. 13 Comparison of specific heat rejection of parallel and 
serial radiators with variation in tube diameter 


comparison between the two configurations for optimum geometric 
and flow parameters for given heat rejection is addressed in 
Sec. 4.1. 


4.1 Optimization and Comparison of Optimum Radiator 
Configurations. The Taguchi-Neuro-Genetic approach was 
applied for both parallel and serial radiator configurations. The 
objective function chosen is the composite objective function, 
which is the sum of nondimensional mass and nondimensional 
pumping power with weighing factors as given in Eq. (8), with 
equality constraint imposed on the total heat rejection. In addition 
to the equality constraint on total heat rejection, inequality con- 
straints on the Pitch of the tubes based on the fabrication aspects 
were also considered. The minimum bend radius of Aluminum 
tubes shall be greater than three times the outer diameter of the 
tube [24]. Considering this aspect, a constraint is posed on the 
pitch of tubes in serial radiator configuration. In parallel radiator 
configuration, the tubes have to be welded to the inlet and outlet 
headers, for which a minimum clearance of 30 mm between the 
tubes is required. Accordingly, the pitch of the tubes in parallel radi- 
ator configuration are constrained. 

Three different cases were considered for optimization study. For 
case 1, equal weightage is given for both mass and pumping power 
by choosing w, and wy, as |. For case 2, full weightage is given to 
mass by choosing w; = 1 and w2=0. For case 3, full weightage is 
given to pumping power by choosing w; =0 and w2=1. 

Tables 7 and 8 provide the input parameters and the results of 
optimization, respectively. For case 1, where equal weightage is 
given to minimize both mass and pumping power, the optimum 
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Fig. 14 Comparison of performance of parallel and serials radi- 
ators with variation in fin thickness 
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Fig. 15 Comparison of specific heat rejection of parallel and 
serials radiators with variation in fin thickness 
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P-Q: Heat rejected by the parallel radiator; 
S-Q: Heat rejected by the serial radiator; 


P-M: Mass of the parallel radiator; 
S-M: Mass of the serial radiator; 


Fig. 16 Comparison of performance of parallel and serial radia- 
tors for varying pitch of tubes 


parallel radiator is about 20% heavier and requires about 35% more 
pumping power than the optimum serial radiator. As the objective 
of the optimization is to minimize the composite objective function, 
the value of composite objective function, which is about 27% 
higher for parallel radiator than the serial radiator, indicates that 
the serial radiator is a better performer. 

For case 2, where total weightage is given to minimize the mass 
of the radiator, the optimum parallel radiator is about 102% heavier 
than the optimum serial radiator. For case 3, where total weightage 
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Fig. 17 Comparison of specific heat rejection of parallel and 
serial radiators with variation in pitch of tubes 
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P-II: Pumping power for the parallel radiator; 
S-II: Pumping power for the serial radiator; 


P-Q: Heat rejected by the parallel radiator; 
S-Q: Heat rejected by the serial radiator; 


Fig. 18 Comparison of performance of the parallel and serial 
radiators for varying pitch of tubes 


Table 7 Input parameters for optimization 
Total heat rejected 270 W 
Area of radiator 1.5x1m 
Fluid inlet temperature, 7;, 283 K 


is given to minimize the pumping power requirement of the radia- 
tor, the optimum parallel radiator requires about 64% higher 
pumping power than the optimum serial radiator. 

For all three optimization cases, the serial radiator displays supe- 
rior performance in terms of lower mass and lower pumping power 
requirement for a given heat rejection than the parallel radiator. 
From the performance analysis, it was observed that at any mass 
flowrate, serial radiator requires higher pumping power than the 
parallel radiator. However, from the optimization study, it is clear 
that for a given heat rejection, mass flowrate requirement of serial 
radiator is lesser than the parallel radiator, causing the serial radiator 
to consume less pumping power than the parallel radiator. 

The optimum configurations obtained from the Taguchi-Neuro- 
Genetic approach are analyzed through conjugate heat transfer 
simulations to verify the output. It was observed that there is a 
good match between the values of heat rejection, mass, and 
pumping power requirement obtained from the conjugate heat trans- 
fer simulations and those predicted by ANN. 


Table 8 Optimization results 


Case 1: wy=1; Case 2: w,=1; Case 3: w, =0; 
W2> 1 W2>= 0 W2> 1 
Parallel Serial Parallel Serial Parallel Serial 
Inner diameters (d;), 4 8.7 4 4.0 4.0 8.7 
mm 
Fin thickness (f)), 1.2 0.6 1.4 0.6 1.6 1.6 
mm 
Pitch of tubes (P;), 102.5 108.0 135 180.0 30.0 108.0 
mm 
Mass flowrate in 50 32 100 70.0 55.0 27 
radiator (Mot), g/S 
Total mass of Tale 6.45 7.27 3.59 14.36 10.11 
radiator (M), kg 
Pumping power (II), 7.64 5.67 17.97 340.30 7.34 4.47 
W 
Value of objective 15.36 12.12 7.27 3.59 7.34 4.47 
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Fig. 19 Temperature palettes of optimum radiator configurations: (a) parallel radiator and 


(b) serial radiator 


Figures 19(a) and 19(b) display the temperature palette obtained 
from conjugate heat transfer simulations for optimum parallel and 
serial radiator configurations, respectively. The hottest zone in the 
parallel radiator is at the top of the radiator, and the coldest zone 
is at the bottom of the radiator. The significant temperature gradient 
is along the tube length, from top to bottom of the radiator, com- 
pared with the gradient along the fin length. In contrast, the tem- 
perature gradient is along the width of the radiator in serial 
radiator configuration because of the fact that the flow enters at 
the left edge, travels all along the radiator flow path and exits at 
the right edge. 


5 Conclusions 


A serial radiator configuration, which has proven advantages in 
the ground applications, is proposed and analyzed for space appli- 
cations. Parallel radiator configuration, which is well studied in 
the literature, is taken as a baseline configuration for comparing 
the performance of serial radiator configuration. Conjugate heat 
transfer simulations were performed by varying one variable param- 
eter at a time to carry out performance analysis. From the perfor- 
mance analysis, it was observed that the specific heat rejection 
capacity, which is a ratio of total heat rejection to the mass of radi- 
ator, is higher for the proposed serial radiator than the conventional 
parallel radiator for any tube diameter, the thickness of fin and pitch 
of the tubes under consideration. The pumping power required by 
the serial radiator is higher than the parallel radiator for any mass 
flowrate under consideration. However, the total heat rejection is 
lower for parallel radiator than serial radiator for any mass flowrate. 
Thus, optimization of radiators is carried out by constraining the 
total heat rejection. 

The Taguchi-Neuro-Genetic approach is adopted for performing 
design optimization. The conjugate heat transfer simulations were 
carried out for the partial factorial set of simulations determined 
by Taguchi’s orthogonal arrays, covering the entire design 
domain. The simulation results were used to train ANN. The 
trained ANN predicts the output with an error of less than 1%. 
The trained ANN was coupled with the genetic algorithm to carry 
out multi-parametric optimization of radiator configurations. 
Through the Taguchi-Neuro-Genetic approach, the optimum paral- 
lel and serial radiator configurations are obtained. Optimization is 
conducted for three different cases. For case 1, where equal weigh- 
tage is given to minimize both mass and pumping power for given 
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heat rejection, the optimum parallel radiator is about 20% heavier 
and requires about 35% more pumping power than the optimum 
serial radiator. For case 2, where total weightage is given to mini- 
mize the mass of the radiator for given heat rejection, the 
optimum parallel radiator is about 102% heavier than the 
optimum serial radiator. For case 3, where total weightage is 
given to minimize the pumping power requirement of the radiator 
for given heat rejection, the optimum parallel radiator requires 
about 64% higher pumping power than the optimum serial radiator. 
From the results of case 1 and case 3, when an equality constraint is 
posed on total heat rejection, the pumping power requirement of the 
serial radiator is lower than the parallel radiator because the serial 
radiator requires a lower mass flowrate than the parallel radiator 
to reject a given heat load. From all three optimization cases, the 
proposed serial radiator turns out to be a superior performer 
in terms of lower mass and lower pumping power requirement 
than the conventional parallel radiator. Thus, it is concluded that 
for the spaceflights, where mass and power are of prime importance, 
the proposed serial radiator configuration has proven advantageous 
over conventional parallel radiator configuration. The serial radiator 
configuration and the optimization study results presented in this 
study can be applied to radiator design for human space flights 
and high-power satellites. 
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Nomenclature 


h = heat transfer coefficient 

n = number of variables 

A = area of radiator exposed surface 

F = view factor from the surface to deep space 

L = length of the radiator 
M = total mass of the radiator including metal and fluid 

masses per unit area 
Q = heat rejected by the radiator in Watts 
E 


= Earth Shine 
d; = inner diameter of the tube in the radiator 
d;, = diameter of tube used for nondimensionalization 
l,; = level of variance of ith variable 
Mdot = Mass flowrate of coolant in radiator 
Mdot,o = Mass flowrate used for nondimensionalization 
m, = mass of the radiator panel 
mz = mass of tubes in the radiator 
m3 = mass of fluid in the tubes 
m4 = mass of header tube (only for parallel radiator 
configuration) 
ms = mass of the fluid in header tubes (only for parallel 
configuration) 
ext = external heat load on radiator (W/m7) 
te = thickness of radiator fin 
to = thickness of fin used for nondimensionalization 
Ww, = weighing factor for mass 
W2 = weighing factor for pumping power 
C, = specific heat of the coolant/fin material 
M, = mass of radiator used for nondimensionalization, | kg 
Ntaguchi = Minimum number of simulations to be conducted based 
on Taguchi method 
P; = pitch of the tubes in radiator 
Pio = pitch of the tubes used for nondimensionalization 
Ty = average fluid temperature in the tube 
T;, = temperature of coolant at the inlet of the tube 
Taq = average surface temperature of radiator 
Tsink = Sink temperature 
Tspace = deep space temperature 
a, = solar absorptivity of the exposed surface of radiator 
AP = pressure drop in the entire ERL 
¢ = IR emissivity of the radiating surface 
1 = efficiency of the radiator 
Np = efficiency of the pump 
A = Albedo 
Il or PP = pumping power 
Il, = pumping power used for nondimensionalization, | W 
py = density of fluid 
Py = density of radiator tube and fin material 
o = Stefan Boltzmann constant 
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